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Planning  a  Purely  Translational  Motion  for  a  Convex  Object 
in  Two-Dimensional  Space  Using  Generalized  Voronoi 

Diagrams 

Daniel  Leven  and  Micha  Sharir 

i 

School  of  Mathematical  Sciences, 

Tel  Aviv  University 
ABSTRACT 

An  0{n  log  n)  algorithm  for  planning  a  purely  translational 
motion  for  a  simple  convex  object  amidst  polygonal  barriers  in 
two-dimensional  space  is  given.  The  algorithm  is  based  on  a  new 
generalization  of  Voronoi  diagrams  (similar  to  that  proposed  by 
Chew  and  Drysdale  [CD]  and  by  Fortune  [Fo]),  and  adapts  and 
uses  a  recent  technique  of  Yap  for  the  efficient  construction  of 
these  diagrams. 

1.   Introduction 

In  this  paper  we  consider  the  following  problem.  Given  a  convex  object 
B  of  relatively  simple  shape  and  a  set  of  polygonal  obstacles,  which  are 
composed  of  a  set  of  n  line  segments  ("walls"),  plan  a  purely  translational 
motion  of  B  from  a  given  initial  position  to  a  desired  final  position  during 
which  B  avoids  collision  with  the  obstacles,  or  determine  that  no  such  motion 
exists.  A  recent  paper  by  Kedcm  and  Sharir  [KS]  (see  also  [KLPS])  has 
presented  an  0{n  log-  /;)  time  algorithm  for  solving  this  problem,  based  on 
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calculating  the  vector  differences  At  —  B  of  each  of  the  obstacles  A,  and  B, 
and  then  on  calculating  the  union  of  these  differences,  obtaining  an  "obstacle 
space"  K  for  a  reference  point  on  B,  and  thus  reducing  the  motion-planning 
problem  to  the  calculation  of  the  connected  components  of  the  complement  of 
K.  Another  related  work  is  that  of  O'Dunlaing  and  Yap  [OY]  for  the  specific 
case  where  5  is  a  disc.  They  showed  that  if  the  free  space  available  for  B  is 
bounded,  a  motion  of  B  between  two  positions  Z,  Z'  exists  if  and  only  if 
there  exists  a  motion  of  (the  center  of)  B  along  the  Voronoi  diagram  of  the 
obstacles  from  some  initial  position  W  to  some  final  position  W  on  that 
Voronoi  diagram,  where  W,  W  are  suitably  defined  retractions  of  Z,  Z'  into 
the  Voronoi  diagram.  The  complexity  of  the  algorithm  in  [OY]  is  0{n  log  n), 
which  is  the  time  required  to  construct  the  Voronoi  diagram  for  a  set  of  n 
line  segments  [Ya]. 

The  idea  behind  the  present  paper  is  similar  to  that  of  [OY].  The  reason 
that  the  Voronoi  diagram  is  a  suitable  tool  for  the  case  where  5  is  a  disc  is 
that  given  a  position  Z  of  (the  center  of)  B  which  is  not  on  the  Voronoi 
diagram  of  the  obstacles,  we  can  move  B  in  a  canonical  manner  to  a  position 
W  on  the  Voronoi  diagram  in  such  a  way  that  the  the  distance  of  B  from  the 
set  of  obstacles  when  positioned  at  W  is  greater  than  the  distance  of  B  from 
the  obstacles  when  positioned  at  Z.  Hence  no  harm  is  done  if  we  move  B 
from  Z  to  W.  Our  idea  is  to  define  a  generalized  Voronoi  diagram  which  has 
a  similar  property  for  a  given  convex  object  B  of  relatively  simple  shape  (e.g. 
a  polygon  with  a  small  number  of  sides)  which,  however,  is  allowed  only  to 
translate  and  not  to  rotate  (for  a  disc  this  further  constraint  is  of  course 


redundant).  As  it  turns  out,  this  generalized  diagram  has  similar  properties  to 

those  of  more  conventional  Voronoi  diagrams;  it  has  0(n  +  A^)  edges  and 

vertices,  and  can  be  constructed  using  similar  techniques  in  0{n  log  N)  time, 

where  N  is  the  number  of  obstacles  and  n  is  the  total  number  of  obstacle 

corners,  and  where  the  complexity  of  B  is  assumed  to  be  fixed.  This  in  turn 

also  yields  an  0(n  log  N)  time  algorithm  for  planning  a  purely  translational 

i 
motion  of  arbitrary  convex  objects  of  simple  shape. 

The  idea  of  using  generalized  Voronoi  diagrams  of  the  sort  proposed 
here  has  recently  and  independently  been  suggested  by  Chew  and  Drysdale 
[CD]  and  by  Fortune  [Fo].  However  they  had  different  motivations  for 
introducing  such  a  structure;  for  example,  Chew  and  Drysdale  were 
interested  in  the  problem  of  determining  the  largest  convex  shape  homothctic 
to  a  given  shape,  which  will  fit  in  the  plane  without  touching  any  of  n  given 
points.  Their  analysis  thus  corresponds  to  assuming  that  the  given  obstacles 
are  all  singleton  points.  The  present  paper  contains  a  detailed  analysis  of  the 
resulting  diagram,  which  is  more  general  than  the  one  studied  in  [CD], 
presents  a  different  method  for  the  calculation  of  the  diagram,  and  discusses 
applications  of  the  diagram  to  the  translational  motion  planning  problem 
mentioned  above. 

This  paper  is  organized  as  follows.  In  section  2  we  define  the  appropriate 
generalization  of  Voronoi  diagrams  which  we  call  fi-Voronoi  diagram,  and 
state  and  prove  its  main  properties.  In  section  3  we  present  an  algorithm  for 
the  calculation  of  this  diagram  in  0(n  log  A^)  time  and  in  Section  4  we  show 
how  to  use  this  diagram  for  the  solution  of  the  above  mentioned  translational 
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motion  planning  problem. 

2.  The  B-Voronoi  Diagram  for  a  Set  of  Convex  Polygonal  Objects 

Let  5  be  a  compact  convex  object  of  some  simple  shape  (the  term  simple 
will  be  defined  more  precisely  in  the  next  section),  which  we  assume  to  have 
non-empty  interior,  and  let  O  be  an  arbitrarily  chosen  fixed  interior  point  of 
fi.  Let  5  be  a  set  of  A^  closed  interior-disjoint  polygonal  objects  ("obstacles") 
composed  of  a  total  of  n  line  segments  ("walls")  which  form  the  boundary  of 
the  free  space  V  available  for  B.  We  assume  that  each  object  in  S  is  convex 
and  allow  objects  in  S  to  share  a  common  wall  (in  this  way  non  convex 
obstacles  may  be  represented  as  the  union  of  several  interior-disjoint  convex 
subparts).  To  avoid  certain  technical  difficulties,  we  will  assume  that  the 
obstacles  in  S  and  B  arc  placed  in  "general  position"  and  thus  exclude 
degenerate  configurations  of  these  objects.  These  assumptions  are  required 
either  to  prove  the  properties  which  we  state  in  this  section  or  for  the 
algorithm  which  we  describe  in  section  3,  and  are  as  follows. 

Assumptions  on  General  Position. 

(a)  No  boundary  edge  of  i?  or  of  an  obstacle  in  S  is  vertical. 

(b)  No  boundary  edge  of  B  is  parallel  to  the  boundary  edge  of  an  obstacle  or 
to  a  line  joining  a  pair  of  boundary  corners  of  these  obstacles. 

(c)  No  more  than  three  objects  in  S  can  be  tangent  to  any  homothetic  copy  of 
B  and  lie  in  its  exterior. 

(d)  No  homothetic  copy  of  B  is  such   that  two  adjacent  boundary  segments  of 


it  simultaneously  touch  two  wall  corners  and  their  common  endpoint  touches 
a  wall. 

Remarks:  (1)  Assumption  (a)  is  required  only  for  the  construction  of  the 
diagram  in  Section  3,  but  not  for  proving  its  properties,  and  can  always  be 
ensured  by  an  appropriate  rotation  of  the  coordinate  frame. 

(2)  Assumption  (b)  is  required  to  ensure  that  the  5-Voronoi  diagram  will  be 
one-dimensional.  Degenerate  configurations  which  occur  when  this 
assumption  is  not  satisfied  can  be  handled  by  standard  techniques;  for 
example  one  can  use  the  method  of  infinitesimal  perturbations  mentioned  in 
[SS]. 

(3)  Assumption  (c)  resembles  the  assumption,  usually  made  in  the  analysis  of 
the  standard  Voronoi  diagram  for  a  set  of  points,  that  no  more  than  three 
points  in  the  set  may  be  cocircular  (cf.  [Sh]).  If  this  assumption  is  violated 
the  diagram  may  have  vertices  of  degree  greater  than  three.  This  will  cause 
some  problems  in  the  contour  tracing  procedure  described  in  section  3  but,  as 
in  the  calculation  of  standard  diagrams,  the  algorithm  can  be  slightly 
modified  to  handle  such  a  degeneracy. 

(4)  Assumption  (d)  is  required  to  prove  property  (9)  below  which  states 
roughly  that  the  diagram  does  not  contain  "degenerate"  vertices.  In  practice  a 
slight  modification  of  our  algorithm  will  suffice  to  make  it  work  even  if  this 
assumption  is  not  satisfied. 

We  will  refer  to  the  position  of  B  in  the  plane  in  which  O  lies  at  the 
origin  as  the  standard  position  of  B  and  denote  by  Bq  the  set  of  points  in  B 
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when  given  that  standard  position.  Let  p ,  q  he  two  points  in  V.  We  define 
the  B -distance  from  p  Xo  q  (or  equivalently,  the  B -closeness  of  p  to  ^), 
denoted  hy  dgip ,  q),  os. 

Note  that  since  O  is  an  interior  point  of  B,  the  5-distance  function  is  always 
finite  and  continuous,  obeys  the  triangle  inequality  but  need  not  necessarily 
be  symmetric,  and  is  thus  not  induced  by  a  metric  in  general  (such  a  notion 
was  introduced  by  Minkowski  in  1911  and  was  called  a  convex  distance 
function  (cf.  [La])).  The  5 -distance  from  a  point  p  €  V  to  a  given  object 
Si  i  S  is  defined  in  a  natural  way  as 

ds(p,si)  =  inf{K  :p  +  kBo  D  s^  ¥^  0} 
(Note  that  under  our  assumptions  that  B  and  S  are  in  general  position,  there 

exists  a  unique  y  ^  s,  al  which  the  5-distance  from  p  to  si  is  attained). 

Using  dg  we  define  the  B-Voronoi  diagram  VorB{S)  associated  with  the 
set  S  and  B  as  follows.    For  each  ii^j  define 

H(iJ)^{y  ^  £2  .  ^^^^^^)  <  ds(y,sj)}, 
i.e.  the  set  of  all  points  whose  5-distance  to  si  is  no  greater  than  their  B- 

distance  to  sj.    Then  define  the  (closed)  B-Voronoi  cell  Cs{si)  associated  with 

Sj  to  be 

i.e.  the  set  of  all  points  y  whose  5-closeness  to  Si  is  no  greater  than  v's  B- 
closcncss  to  any  other  clement  of  S.  Finally,  the  B-Voronoi  diagram  Vorg{S) 
is  defined  to  be  the  set  of  points  which  belong  to  more  than  one  i?-Voronoi 
cell.    The  edges  of  this  diagram  arc  maximal  connected  arcs  which  belong  to 
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exactly  two  B-Voronoi  cells  and  the  vertices  are  those  points  which  belong  to 
more  than  two  cells.  To  simplify  the  following  analysis  we  shall  assume  that 
for  each  pair  of  objects  s^,  $2  in  S  which  share  a  common  boundary  segment 
/,  the  exterior  angles  formed  between  Si,  si  at  the  endpoints  of  /  are  both  less 
than  or  equal  to  180°,  and  thus  the  locus  of  points  equidistant  from  any  given 
pair  of  objects  in  S  is  always  a  one-dimensional  set  (if  that  angle  is  180°,  this 
property  follows  from  our  assumption  that  B  and  S  are  in  general  position), 
so  that  in  this  case  the  5-Voronoi  diagram  is  indeed  one-dimensional.  If  this 
assumption  is  violated  at  some  point  x  common  to  the  boundary  of  two 
obstacles  s^,  si  (i.e.  the  exterior  angle  between  si  and  si  at  x  is  greater  than 
180°),  then  we  can  treat  j:  as  a  singleton  obstacle,  replace  s-^,  Si  by  Si  —  [x] 
,S2  —  {x}  respectively,  and  redefine  5-closeness  similar  to  the  definition  of 
closeness  in  [Ki],[Ya].  However,  for  the  sake  of  simplicity  we  will  continue 
to  make  this  assumption,  which  can  always  be  enforced  when  partitioning  a 
non-convex  obstacle  into  a  collection  of  openly-disjoint  convex  parts.  By  our 
assumptions  that  B  and  S  are  in  general  position  no  point  on  Vorfi(5)  lies  in 
more  than  three  5-Voronoi  cells.  Fig.  2.1  shows  an  example  of  such  a  B- 
Voronoi  diagram.  Note  that  the  5-Voronoi  diagram  for  the  case  where  Z?  is  a 
disc  is  exactly  the  standard  Voronoi  diagram  of  the  set  S  of  obstacles  [OY]. 

The  5-Voronoi  diagram  just  defined  has  the  following  properties. 

(1)  The  collection  of  5-Voronoi  cells  covers  the  whole  plane. 

Proof:  Follows  immediately  from  the  definition  of  d^  and  the  fact  that  O  is 
an  interior  point  of  B. 
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Fig.  2.1  TheB-Voronoi  diagram. 

(2)  Let  J  be  a  point  in  some  cell  Cb(5,)  and  let  y'  be  the  B-closest  point  to  it 
on  Si-  Then 

(a)  The  line  segment  I  ^  yy'  is  wholly  contained  in  Cb{s,)  and  in  addition,  if 
there  is  a  point  z  £  /  which  is  contained  also  in  another  cell  Cg{sj)  then  the 
whole  segment  yz  is  also  contained  in  C^isj)  (we  capture  this  property  by 
calling  Cg^Si)  star  shaped  with  respect  to  Sj.) 

(b)  If  p  is  any  point  such  that  y'yp  is  a  straight  line,  then  the  B-closest  point 
on  Si  to  p  is  y' . 

Proof:  (a)  Suppose  that  there  is  an  interior  point  z  C  /  which  belongs  to  a 
second  cell  C^^Sj)  and  let  z'  be  the  5-closcst  point  to  z  on  Sj  (see  Fig  2.2(a)). 
We  have  by  our  assumption 

(^h{^,Sj)  <  ds{z,Si) 


XBo  +  y 


(a) 


Fig  2.2 


(b) 


Let  dB(y,Si)  —  X.  If  the  inequality  were  strict  (i.e.  z  did  not  belong  to 
Cb{si))  then  z'  would  belong  to  the  interior  of  the  triangle  yy'y"  where  yy" 
is  a  line  segment  parallel  to  zz'  such  that  dgiyj")  —  X.  But  since  the 
triangle  yy'y"  is  wholly  contained  in  y  +  XBQ,  this  would  imply 
^aiyj^j)  ^  ^siyy^i)  contrary  to  our  assumptions.  Therefore  z'  must  lie  on 
the  segment  y'y"  and  also  be  a  boundary  point  of  y  +  \BQ  which  however 
implies  that  every  point  on  /  between  y  and  z  is  equally  close  to  Si,  Sj  as 
asserted. 

(b)  Let  X  be  the  scalar  multiple  of  B  such  that  XBq  +  y  touches  si  at  v' .  Since 
y'  is  the  closest  point  to  y  on  Si  and  since  B  and  Si  are  convex  there  is  a  Hne  T 
passing  through  y'  which  separates  them  (see  Fig.  2.2(b)).  Let  X'  >  X  be 
the  multiple  of  B  such  that  X'Bq  +  p  also  touches  Sj  at  y' .  It  is  easily 
checked  that  T  also  separates  X'Bq  +  p  and  5,.  But  then  since  B  and  Sj  arc 
convex  it  follows  that  the  fi-distance  from  p  to  any  point  on  T  (and  hence 
also  to  any  point  in  Si)  is  greater  than  or  equal  to  X',  and  this  proves  our 
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claim. 

(3)  Suppose  1 1,  I2  are  two  line  segments  in  Cb{si)  which  connect  two  points 
xi,  x-i  in  Cfl(5'/)  to  their  respective  B-closest  points  y^,  y^  on  5,.  Then  either 
1 1, 1 2  do  not  intersect  except  possibly  at  their  endpoints  or  else  one  of  these 
segments  is  wholly  contained  in  the  other. 

Proof  (The  technique  used  here  resembles  that  of  [KS]):  Suppose  /j,  1 2 
intersect  at  a  single  point  and  let  dB{xi,  Si)  —  \,  dg{x2,  St)  =  X'.  Let  z^ 
(resp.  Z2)  be  a  point  at  a  B-distance  X  (resp.  X')  from  xi  (resp.  X2)  in  the 
direction  X2y2  (resp.  x-^yi)  (see  Fig.  2.3). 


Clearly  the  triangle  x^z-^i  is  contained  in  x^+XEq  and  the  triangle  X2Z2V2  is 
contained  in  X2  +  >^'Bq  and  therefore  both  these  triangles  do  not  contain  a 
point  (in  particular  the  points  y^,  y2)  of  Sj  in  their  interior.  But  since  the  two 
triangles  aiz^Vi,  X2Z2y2  are  similar  and  have  pairwise  parallel  sides,  and  their 
two  sides  /j,  I2  intersect  one  another,  we  must  have  either  y^  ^  Aj:2Z2.V2  or 
V2  ^  Aa;iZlVj.  ITicrefore  Zi>'2>'iZ2  must  be  a  straight  line  segment  and  yi>'2 
must  therefore  be  a  subsegmcnt  of  the  boundary  of  5/.  But  this  means  that 
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the boundary  of  B  contains  a  line  segment  which  is  parallel  to  a  boundary 
segment  of  j,  which  contradicts  our  assumptions  of  general  position. 

(4)  Each  edge  of  the  fi-Voronoi  diagram  is  simple  (i.e.  the  removal  of  a 
single  point  on  such  an  edge  divides  it  into  at  most  two  connected 
components). 

Proof:  Let  p  be  a  point  on  an  edge  e  on  the  boundary  of  Cg^si),  €^{52)  such 
that  the  removal  of  p  divides  e  into  more  than  two  disconnected  components, 
and  let  C^,  Ci,  C3  be  any  three  such  components  in  some  small  neighborhood 
of  p.  Define  4>/(a:),  /  =  1,2,  for  j:  ^  e  to  be  the  5-closest  point  on  j,  to  x. 
Denote  <^i{p)  by  pi.  Plainly  <^i  is  continuous  and  thus  ^i(Ci^),  k  =  1,2,3,  are 
connected  portions  on  the  boundary  of  s^  whose  closure  contains  p-^.  llie 
star-shapedncss  of  Cs{si)  and  the  fact  that  the  5-diagram  is  one-dimensional 
easily  imply  that  pi  is  an  endpoint  of  all  three  projections  ^liC),)  and  that 
these  projections  are  all  disjoint.  Thus  one  of  them,  say  ct)i(Ci),  must  be  the 
singleton  point  p^.  But  then  property  (2b)  implies  that  C^  must  be  a  line 
segment  all  of  whose  points  are  equally  5-closest  to  p^  and  to  S2.  Arguing  as 
in  the  proof  of  (2),  it  follows  that  p-jJi  must  be  parallel  to  a  side  of  B;  thus 
our  assumptions  on  general  position  of  B  and  S  imply  that  p2  must  be  a 
corner  of  5^2  (see  Fig.  2.4).  However,  a  symmetrical  argument,  using  ^2 
instead  of  ct)i  would  imply  that  p-^  must  be  a  corner  of  s^,  thereby 
contradicting  our  assumptions  of  general  position,  namely  that  no  side  of  B 
is  parallel  to  the  line  joining  two  corners  of  obstacles  in  S. 

(5)  Let  a  be  a  point  on  the  locas  C  of  points  which  are  equally  5-close  to  two 
distinct  objects  5,,  Sj  and  let  >•,■  (rcsp.  yj)  be  the  (unique)  5-closest  point  to  a 
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c^^  (C2) 


Fig  2.4 

on  si  (resp.  Sj).  Denote  the  fi -distance  from  a  to  yi  and  yj  hy  \.  Let  7,  (resp. 
Tj)  be  a  common  tangent  at  yi  (resp.  y^)  to  5/  (resp.  Sj)  and  c  +  XBq  and 
suppose  that  these  tangents  are  unique.  Then  the  (unique)  tangent  to  C  at  a 
is  the  line  T  which  passes  through  a  and  the  intersection  point  z  of  Tt  and  Tj. 

Proof:  Let  a'  be  a  point  on  T  sufficiently  near  a,  let  y/  be  a  point  on  Ti  such 
that  a'yi'  is  parallel  to  ay^,  and  let  yj'  be  a  point  on  Tj  defined  in  an 
analogous  manner.   Let 


X'  =  X- 


P'yi 


Then  it  is  easily  seen  that  a'  +  X'Bq  contains  >•;',  y/  and  is  fully  contained  in 
the  wedge  between  T,  and  Tj.  Since  at  least  one  of  the  two  objects  Si  (resp. 
Sj)  and  a  +  XBq  is  smooth  at  their  point  of  contact  v,-  (resp.  yj),  it  follows 
that  the  5-distances  of  a'  from  both  Sj  and  Sj  are  equal  to  X'  up  to  first  order 
approximation,  proving  that  T  is  indeed  the  tangent  to  C  at  a.  This  argument 
can  also  be  used  to  show  that  T  is  the  unique  such  tangent. 
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(6)  Corollary:  (a)  If  the  boundary  of  B  is  smooth  (i.e.  B  has  a  unique 
tangent  at  each  point  on  its  boundary)  then  each  5-Voronoi  edge  is  smooth. 

(b)  If  5  is  a  polygonal  object  then  each  edge  of  the  5-Voronoi  diagram  is  a 
polygonal  arc.  Furthermore,  if  B  has  k  sides  then  the  total  number  of  comers 
in  the  edges  of  Vorg{S)  is  at  most  kN  +  n. 

Proof:  (a)  is  trivial,  because  in  this  case  the  conditions  of  (5)  hold  at  each 
interior  point  of  any  of  the  edges  of  the  diagram. 

(b)  Let  e  be  a  5-Voronoi  edge  separating  Cb(si)  from  Cg(sj).  Clearly  e  can 
be  broken  into  portions  at  each  of  which  the  two  contacts  between  the 
appropriate  expansion  of  B  and  s,,  Sj  are  made  between  fixed  corners  or 
segments  of  B  and  fixed  segments  or  corners  of  s,  and  sj.  As  long  as  these 
type  of  contact  are  fixed,  it  is  easily  seen  that  e  traces  a  straight  line  segment 
(this  follows  e.g.  from  (5)).  Moreover,  this  straight  segment  is  terminated 
when  contact  is  made  between  a  corner  of  B  and  a  corner  of  either  Si  or  sj.  It 
can  be  shown  that  if  O  moves  around  the  boundary  of  Csisi)  with  B 
expanding  or  shrinking  so  as  to  maintain  contact  with  Si  then  the  total 
number  of  points  at  which  a  corner  of  B  makes  contact  with  a  comer  of  Si 
during  that  traversal  is  at  most  ^  +  n,,  where  n^  is  the  number  of  sides  of  Si 
(cf.  e.g  [GRS]).  Hence,  summing  over  all  5-cells  of  the  diagram,  we  obtain 
the  asserted  bound. 

Remark:  The  same  arguments  used  to  prove  (6b),  combined  with  the 
observations  in  [GRS],  can  be  used  to  show  that,  under  the  assumptions  in 
(6b),  the  locus  of  points  5-equidistant  to  two  objects  Sj,  Sj  having  /i,,  rij  sides 
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respectively,  can  be  calculated  in  time  0{k  +  n,  +  rij).  This  observation  will 
be  significant  in  the  algorithm  for  constructing  the  diagram,  given  in  section 
3  below. 

(7)  VorB{S)  need  not  be  connected.  In  fact  it  can  have  up  to  0{n)  connected 
components.  However,  every  connected  component  of  VorB{S)  is  unbounded. 

Proof:  Since  the  5-Voronoi  diagram  for  B  a  disc  is  the  standard  Voronoi 
diagram  of  5",  the  first  claim  follows  from  the  fact  that  the  standard  Voronoi 
diagram  may  have  0{n)  connected  components  (cf.  [Sha]).  The  diagram  can 
also  become  unbounded  due  to  nonsmooth  comers  of  B.  Fig.  2.5  shows 
examples  of  both  these  types  of  unboundedness. 


Fig  2.5 

Let  AT  be  a  bounded  connected  component  of  VoreiS).  Since  the  set  of  B- 
closeness  cells  covers  the  whole  plane  there  is  an  object  si  whose  cell 
surrounds  K.  Obviously  Si  lies  outside  K.  Let  y  be  any  point  in  K  and  let  x 
be  the  5-closest  point  to  it  on  si.  Continue  the  segment  xy  past  >'  and  away 
from  Si  until  is  reaches  an  interior  point  p  in  C^isi)  (such  a  point  must  exist 
since  Cgisi)  surrounds  A.').    But  then  by  property  (2b)  the  closest  point  to/;  is 
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X  and  this  contradicts  the  fact  that  Cb{si)  is  star-shaped. 

(8)  Corollary:  There  are  no  isolated  points  in  Vors{S). 

(9)  Let  /?  be  a  vertex  of  VoreiS)  on  the  boundary  of  three  cells 
Cb{si),  Cb{sj),  Cgisk)-  Then  every  neighborhood  of  p  intersects  the  interiors 
of  all  three  cells. 

Proof:  Suppose  to  the  contrary  that  there  is  a  point  p ,  whose  fi-distances  to 
Si,  Sj,  Si-  are  all  equal  to  some  X,  such  that  there  exists  a  small  neighborhood 
N  of  p  which  contains  no  interior  point  of  Cb(si).  We  may  assume  that  N 
contains  no  point  of  a  fourth  cell,  since  if  this  were  the  case  for  all 
sufficiently  small  such  neighborhoods  N,  then  p  would  be  equally  B-close  to 
four  objects  and  we  have  ruled  out  such  cases  by  our  assumptions  on  general 
position.  Let  the  B-closest  point  to  p  on  Sj  be  >>,■  and  consider  the  intersection 
/  of  the  line  segment  /?>',  with  A^  By  property  (2)  pyi  is  contained  in  ^(5',) 
and  therefore  it  cannot  intersect  interior  points  of  any  other  cell.  Also,  by 
our  assumption,  there  are  no  interior  points  of  Cs(.s'/)  on  either  side  of  a 
portion  of  /  near  p ,  and  therefore  that  portion  of  /  lies  entirely  in  the  interior 
of  the  union  of  Cb{sj)  and  Cb{s)^).  But  (as  in  the  proof  of  property  (2))  since 
B  and  the  objects  in  S  are  all  convex,  this  is  possible  only  if  yjyi  and  yj^yi  both 
lie  on  the  boundary  of  p  +  \Bq  where  yj  (resp.  >'^)  is  the  5-closest  point  to  p 
on  Sj  (resp.  s;^)  (see  Fig.  2.6).  This  however  contradicts  our  assumptions  that 
B  and  S  are  in  general  position  (assumption  (d)  at  the  beginning  of  section 
2),  and  thus  establishes  our  claim. 

(10)  Vor^^S)  consists  of  at  most  0{N)  edges  and  vertices,  where  the  edges  are 
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Fig  2.6 

maximally  connected  components  of  Vors{S)  which  lie  on  the  boundary  of 
the  same  two  cells  of  the  diagram  and  the  vertices  are  points  which  are 
equally  5-closest  to  more  than  two  obstacles. 

Proof;  Property  (9)  implies  that  each  neighborhood  of  a  vertex  v  of  the  B- 
Voronoi  diagram  intersects  three  cells  of  the  diagram  and  thus  contains  at 
least  three  edges  separating  these  cells.  Thus  the  degree  of  v  is  at  least  3 
(therefore  by  our  assumptions  the  degree  of  v  is  exactly  3).  From  this  and 
Euler's  formula  for  planar  graphs  we  obtain  the  required  result. 

Remarks:  (a)  Note  that  in  (10)  we  have  taken  each  edge  of  Vorg(S)  to  be  a 
maximal  connected  component  of  Vorg(S)  which  lies  on  the  boundary  of  the 
same  two  cells  of  the  diagram.  In  practice,  e.g.  when  B  is  polygonal,  one 
may  want  to  break  such  an  edge  into  several  (connected)  subparts,  each 
corresponding  to  different  corners  and  edges  of  the  obstacles  in  S  and  of  B 
which  touch  one  another  at  the  appropriate  expansion  of  B.  It  is  not  hard  to 
sec,  using  corollary  (6),  that  if  the  boundary  of  B  is  regarded  as  a  single 
closed  curve,  then  Vor,j{S)  will  consist  of  0(N  +  n)  such  subedgcs.  If  5  is  a 
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polygonal  object  with  k  sides  and  we  wish  to  distinguish  between  subedges  of 
the  diagram  induced  by  different  sides  or  corners  of  B,  then,  as  in  Corollary 
(6b),  Vors{S)  will  consist  of  0{kN  +  n)  subedges,  which,  in  this  specific 
case,  are  all  straight  line  segments. 

(b)  From  (10)  we  can  easily  deduce  that  the  number  of  positions  of  B  in 

which  it  touches  simultaneously   two  distinct  objects  in  S   and  does  not 

I 
penetrate  into  any  other  object  in  S  is  0{N).  Indeed,  any  such  position  must 

lie  on  an  edge  of  Vor^^S).  On  the  other  hand,  any  edge  e  of  VorB{S)  can 

contain  at  most  two  such  positions,  as  follows  from  the  argument  given  in 

[KS].  Thus  the  properties  of  5-Voronoi  diagrams  provide  a  different  proof  of 

this   property.     (A  more  general  proof  wliich  uses  a  somewhat  different 

technique  has  recently  been  given  in  [LiS],  [KLPS].) 

3.   Efficient  Construction  of  the  B-Voronoi  Diagram. 

In  this  section  we  present  an  efficient  algorithm  for  constructing  Vor^  (5) 
for  simple  convex  objects.  By  simple  objects  we  mean  objects  for  which 
certain  operations,  which  are  used  during  the  construction  of  the  5-Voronoi 
diagram  as  described  below,  can  be  accomplished  in  constant  time.  Such 
operations  are:  (a)  computing  the  5-distance  between  any  two  points  in  V  or 
from  a  given  point  in  V  to  a  straight  line;  (b)  calculating  the  intersection 
points  of  the  locus  of  points  which  are  5-equidistant  to  two  sides  or  corners 
of  objects  in  S  and  another  such  curve  or  a  straight  line  segment.  For  more 
complex  objects  the  actual  complexity  of  such  operations  should  be  taken  into 
account  in  estimating  the  complexity  of  the  algorithm. 
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As  will  be  seen,  the  properties  of  VorgiS)  which  have  been  proved  in  the 
previous  section,  and  which  are  quite  similar  to  the  properties  of  more 
conventional  Voronoi  diagrams,  enable  us  to  construct  VorB{S)  by  much  the 
same  technique  as  that  described  by  Yap  [Ya]  for  the  construction  of 
(standard)  Voronoi  diagrams  for  a  set  of  simple  curves. 

Our  algorithm  proceeds  as  follows.  Let  P  be  the  set  of  the  2N  leftmost 

i 

and  rightmost  points  on  objects  in  S  (recall  that  these  points  are  unique  by 
assumption  (a)  of  general  position)  Following  [Ya]  we  divide  the  plane  by 
vertical  lines  into  2N  stabs  such  that  (the  interior  of)  each  slab  contains 
exactly  one  point  of  P.  The  algorithm  runs  in  stages.  In  each  stage  it 
partitions  the  set  of  slabs  into  disjoint  pairs  of  adjacent  slabs  and  combines 
each  such  pair  into  a  new,  larger,  slab.  That  is,  in  the  first  stage  the  2N  slabs 
are  combined  in  disjoint  pairs  into  A^  new  slabs,  in  the  second  stage  these  A^ 

slabs  are  combined  in  pairs  into  —  slabs  and  so  on.  After  0(log  A^)  stages 

only  a  single  slab  remains  and  the  algorithm  terminates. 

With  each  slab  the  algorithm  associates  a  partial  5-Voronoi  diagram  for 
an  appropriate  set  of  objects  which  intersect  that  slab.  These  diagrams  are 
merged  together  whenever  two  adjacent  slabs  are  combined  to  form  a  new 
slab.  However,  to  keep  down  the  cost  of  creating  and  merging  those 
diagrams,  we  will  maintain  for  each  slab  a  partial  diagram  which  covers  only 
those  portions  of  the  slab  that  contain  points  in  P  (see  below  for  a  more 
detailed  description  of  this  process;  note  that  slab  portions  for  which  no 
diagram  is  yet  maintained  arc  simple  "corridors"  extending  through  the  entire 
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width  of  a  slab  and  containing  no  inner  obstacles). 

The  general  merging  step  of  the  algorithm  thus  combines  partial  5- 
Voronoi  diagrams  of  two  adjacent  slabs  into  a  (possibly  partial)  5-Voronoi 
diagram  of  the  combined  slab.  This  is  done  as  follows.  Consider  a  slab  Q 
divided  by  a  vertical  line  into  two  slabs;  a  left  slab  Q^  and  a  right  slab  Qj^. 
The  intersection  of  the  free  spac^  V  with  Q  is  divided  into  several 
components  by  the  boundary  of  objects  which  cross  Q  through  its  entire 
width  from  left  to  right  (i.e.  objects  whose  leftmost  points  lie  to  the  left  of  Q 
and  whose  rightmost  points  lie  to  the  right  of  Q;  see  Fig.  3.1).  Each  such 
component  is  called  a  window;  note  that  our  definition  of  windows  is 
somewhat  similar  to  that  of  quads  in  [Ya]. 

I 

I   I    I 


Fig  3.1 


With  each  window  W  we  associate  a  5-Voronoi  diagram,  defined  over 
the  entire  plane,  which  is  induced  by  the  portions  of  objects  contained  in  W, 
as  follows.    For  each  object  s  which  intersects  (the  interior  of)  W,  let  s"^ 
denote  that  intersection.   In  addition  let   U^y   (rcsp.  D^,)  denote  the  upper 
(resp.  lower)  portion  of  the  slab  Q  bounded  below  (resp.  above)  by  the 
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portion  of  the  object  which  forms  the  upper  (resp.  lower)  boundary  of  W. 
Then  the  fi-Voronoi  diagram  of  W  is  the  fi-Voronoi  diagram  for  the  set  S(W) 
consisting  of  Uw,  Dy/,  and  all  the  nonempty  s^ ,  for  s  ^  S  (see  Fig.  3.2). 


^ 


ca 


w 


'  D 


w 


Fig  3.2 


Note  that  Ut;y,D^,  may  not  be  convex,  in  which  case  they  are  split  into 
convex  subparts  at  their  concave  comers,  each  of  which  is  regarded  as  a 
corner  of  P  (see  the  comment  at  the  beginning  of  section  2).  It  is  clear  that 
the  5-Voronoi  diagrams  of  different  windows  in  the  same  slab  Q  can  overlap 
one  another  only  outside  Q  (in  fact,  as  will  be  noted  below,  such  overlapping 
is  always  of  a  very  trivial  nature),  and  our  algorithm  will  therefore  calculate 
and  maintain  each  of  these  partial  diagrams  separately.  As  already  noted,  we 
will  maintain  5-Voronoi  diagrams  only  for  those  windows  which  contain  at 
least  one  of  the  points  in  P.  This  ensures  that  at  each  stage  of  the  algorithm 
there  are  at  most  2N  windows  for  which  the  diagram  need  be  constructed  (a 
fact  which  is  crucial  for  the  linearity  of  each  stage),  and  also,  of  course,  that 
upon  termination  the  diagram  of  the  last  remaining  slab  will  indeed  be  the 
complete  Z?-Voronoi  diagram  of  S. 


-21- 

The  main  step  that  needs  to  be  described  is  thus  the  step  which  computes 
the  fi-Voronoi  diagram  of  a  given  window  W  of  Q.  W  is  divided  by  the 
vertical  hne  separating  (2l>  Qr  into  two  sets  of  windows,  those  in  the  slab  Qi 
(which  we  call  the  L-windows)  and  those  in  the  slab  Q^  (which  we  call  the 
/^-windows).  Let  L  (resp.  /?)  be  the  set  of  intersections  of  the  objects  in  S 
with  the  L-windows  (resp.  /^-windows).  The  construction  of  the  5-Voronoi 
diagram  of  W  is  accomplished  by  the  following  three  steps.  First  we 
construct  the  partial  diagram  for  each  L-window  (resp.  i?-window)  that 
contains  at  most  one  point  of  P  (and  whose  diagram  has  not  yet  been 
constructed).  Next  we  combine  the  B-Voronoi  diagrams  of  all  the  L-windows 
(resp.  /^-windows)  into  a  full  diagram  VoreiL)  (resp.  VorB{R))  of  the  left 
(resp.  right)  part  of  W.  Finally  we  merge  the  diagrams  of  the  left  and  right 
parts  of  W  into  the  5-Voronoi  diagram  of  W.  We  now  describe  each  of  these 
steps  in  detail. 

Step  1:  Construction  of  the  5-Voronoi  diagram  of  an  L-window  or  an  R- 
window. 

K  the  5-Voronoi  diagram  of  an  L-window  or  an  /^-window  has  not  yet 
been  constructed,  then  either  we  are  in  the  first  (innermost)  stage  of  the 
algorithm  and  the  window  W*  in  question  contains  just  one  point  in  P,  or  W* 
contains  no  point  of  P.  Therefore  at  most  three  objects  s^,  S2,  s^  can  intersect 
W*  in  one  of  the  following  situations;  either  (a)  U^,-  (resp.  D^,..)  consists  of 
(at  most)  two  objects  (it  may  consist  of  two  objects  if  it  is  a  part  of  an 
original   non-convex   object   which   has   been  split   into  convex   subparts   as 
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discussed  at  the  beginning  of  Section  2),  D-^*  (resp.  f/^^,.)  consists  of  one 
object  and  there  is  no  interior  object  in  W* ,  or  (b)  each  of  t/^y.,  D^^.  consists 
of  one  object  and  there  is  one  interior  object  in  W*.  Hence  we  can  compute 
the  B-Voronoi  diagram  of  W*  by  first  computing  the  locus  of  points  which 
are  equally  B-cIose  to  Si,  Sj  (by  Corollary  (6b)  and  the  remark  following  it 
this  can  be  done  in  time  0{mi  +  mj)  if  we  assume  B  to  be  simple  or 
0{k  +  mi  +  m2)  if  we  assume  5  to  be  a  ^-gon,  where  m^  (resp.  m2)  is  the 
number  of  segments  and  points  composing  s^  (resp.  ^'2)),  and  then  by  adding 
to  the  diagram  53.  By  a  careful  implementation,  which  resembles  the  contour 
tracing  procedure  which  is  outlined  below,  this  stage  can  be  done  in  time 
0(mi  -\-  m2  +  m^,)  (or  0(k  +  m^  +  m2  +  m^))  where  m^  is  the  number  of 
segments  composing  the  third  object  s^  (see  also  [US]  for  a  description  of  the 
contour  tracing  technique).  Thus  the  construction  of  the  fi-Voronoi  diagram 
of  W*  requires  0(k  +  m)  time,  where  m  is  the  number  of  comers  of  objects 
in  W*. 

Step    2:    Vertical    merge    of   the    diagrams    of   the   L-windows    (resp.    R- 
windows). 

Let  X,  Y  be  two  adjacent  L-windows  which  are  separated  by  a  portion  of 
an  object  A  and  assume  that  X  lies  above  Y .  We  claim  that  the  only  cells  of 
the  B-diagrams  of  X,  Y  which  can  overlap  are  a  cell  of  the  upper  portion  of  Y 
and  a  cell  of  the  lower  portion  of  X.  This  follows  from  the  following 
argument.  First  observe  that  these  two  diagrams  can  overlap  only  outside  Qi. 
Let  2  be  a  point  lying  (outside  Qi)  in  a  cell  Cti{si)  of  Vorg{L)  where  Si  is  in  X 
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and  is  not  contained  in  Dx-  Then  the  object  5-cIosest  to  z  in  y  must  be  (part 
of)  Uy.  Indeed  suppose  to  the  contrary  that  this  object  5^2  is  not  (part  of)  Uy 
Let  dBiz,si)  =  ki,  dB{z,S2)  =  X2  and  let  p,  ^  St  be  thefi-closest  point  to  z  on 
Si,  i  =  1,2.  Then  it  is  easily  seen  that  X2  ^  Xi  (since  z  C  Cb{si)  in  VorB{L)) 
and  that  both  p^  and  p2  belong  to  z  -i-  X25o.  But  then  the  segment  p-jj2  is 

contained  in  z  +  X2fio  ^^id  it  follows  that  z  +  X25o  must  intersect  the  interior 

i 
of    Uy    (because   Bq   has   no   vertical    boundary   edge),    contradicting    the 

assumption  that  dB{z,S2)  <  dB{z,Uy). 

Thus  to  merge  the  B-diagrams  of  two  adjacent  L-windows  X,  Y  as  above, 
all  we  have  to  do  is  to  intersect  the  cell  of  Dx  in  the  diagram  for  X  with  the 
cell  of  Uy  in  the  diagram  of  Y  to  obtain  the  cell  of  the  object  A  =  Dx  (^  Uy  in 
the  combined  diagram,  and  leave  all  other  cells  in  both  diagrams  unchanged. 
This  can  plainly  be  done  in  time  which  is  linear  in  the  number  of  edges  of  the 
two  cells  Ca(D_x-),  CsiUy)  and  thus  VorB{L)  can  be  constructed  from  the 
partial  diagrams  of  the  L-windows  in  linear  time.  Similarly,  VorBiR)  can  also 
be  constructed  in  a  completely  analogous  manner  from  the  partial  diagrams 
of  the  /? -windows  in  linear  time. 

Step  3:  "Horizontal"  merge  ofVorgiL),  VorB{R). 

The  final  step  merges  VorgiL)  and  VorsiR)  to  obtain  the  5-Voronoi 
diagram  of  5(1^).  Note  that  the  objects  in  S{W)  which  intersect  both  Qi  and 
Q/j  are  split  by  the  vertical  line  separating  (2/,  and  g^  into  two  parts,  one  in  L 
and  the  other  in  R.  The  present  merge  first  calculates  the  5-diagram  for 
L  U  R  (i.e.  with  each  such  object  remaining  split),  and  only  then  patches  up 
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the  diagram  to  obtain  VorB(S(W)). 

To  merge  the  diagrams  VorB{L),  VorgiR),  we  have  to  cx)mpute  the  locus 
C  of  points  which  are  equally  fi-closest  to  an  object  in  L  and  to  an  object  in 
R.  Following  Kirkpatrick  [Ki]  we  call  C  the  contour  separating  L  from  R.  It 
is  easy  to  see  that  C  is  a  simple,  topologically  closed  curve  since  each  vertex 
of  the  combined  diagram  has  degree  3  (for  details  concerning  this  standard 
argument  see  e.g  [LS]).  Hence  C  partitions  the  plane  into  disjoint  open 
connected  regions  each  of  which  consists  either  entirely  of  points  which  are 
5-closer  to  objects  in  L  than  to  objects  in  R  (in  which  case  we  call  such  a 
region,  as  in  [LS],  an  L-region)  or  entirely  of  points  which  are  5-closer  to 
objects  in  R  than  to  objects  in  L  (in  which  case  we  call  such  a  region  an  R- 
region).  Since  by  our  construction  each  object  in  L  lies  to  the  left  of  each 
object  in  R  we  have  the  following  lemma. 

Lemma  3.1:  (Fig.  3.3)  Let  cj,  ,C2,  be  the  common  endpoints  of  segments 

of  objects  in  L  and  in  R  on  the  vertical  line  which  separates  Q^  and  Q/^ ,  and 
suppose  that  these  points  are  given  in  increasing  order  of  their  j -coordinates. 
Then  the  contour  C  consists  of  /  portions  C^,  .  .  .  ,  C,  such  that  C,-  starts  at 
C2,_i  and  ends  at  C2/.  If  the  leftmost  and  rightmost  points  of  B  are  smooth 
then  Cj,  .  .  .  ,  Cj  are  all  simple  connected  arcs;  otherwise  each  of  these 
portions  C,  is  still  simple  but  not  necessarily  connected.  In  this  case  each  Q  is 
the  limit  (in  the  Hausdorff  topology  of  closed  sets)  of  a  sequence  of  simple 
connected  arcs  connecting  the  corresponding  points  C2/-1  and  Co,-;  thus  every 
connected  portion  of  C;  is  unbounded  in  both  directions  (except  at  the 
endpoints  Co,- 1,  c;,). 


Fig.  3.3  The  setting  of  Lemma  3.1 

Proof:  We  will  first  show  that  if  the  the  leftmost  and  rightmost  points  of  B 
are  both  smooth  then  each  pair  of  objects  in  L  can  be  connected  by  an  arc 
which  does  not  intersect  the  contour  and  (by  a  symmetric  argument)  that 
each  pair  of  objects  in  R  can  also  be  connected  by  such  an  arc.  This  will  imply 
that  if  the  leftmost  and  rightmost  points  of  B  are  both  smooth  there  is  only 
one  L-region  and  only  one  /? -region  in  Q  and  hence  the  contour  must  have 
the  form  asserted  by  the  Lemma,  since  it  obviously  passes  through  the  points 
c^,  .  .  .  ,  C2f  For  objects  B  whose  leftmost  or  rightmost  point  is  not  smooth 
we  will  obtain  the  result  by  approximating  B  by  the  vector  sum  B^  of  B  and  a 
disc  with  a  sufficiently  small  radius  e  >  0. 

Details  are  as  follows.  Let  Si,  Sj  €  L  and  let  v,-  (resp.  yj)  be  the  leftmost 
point  on  s,  (resp.  Sj).  For  each  X  >  0  let  a/(X)  (resp.  ay(X))  be  the  leftmost 
point  whose  5-distancc  to  Sj  (resp.  Sj)  is  X.  Since  we  assume  that  B  has  no 
vertical  boundary  edges,  it  follows  that  the  locus  A,  (resp.  Aj)  of  all  points 
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ai{\)  (resp.  aj{\)),  X  ^  0,  is  an  infinite  ray  emanating  leftwards  from  y, 
(resp.  yj),  and  that  Ai  and  Aj  are  parallel.  Let  e^  be  the  line  segment 
connecting  a/(X)  and  aj{k)  and  eo  the  line  segment  connecting  >,  and  yj.  We 
claim  that  if  the  rightmost  point  of  B  is  smooth  then  for  a  large  enough  X 
every  point  p  ^  e^  lies  in  an  L-region,  so  that  St  and  Sj  are  connected  by  the 

polygonal  arc  ytai{k)aj{X)yj  which  is  completely  contained  in  an  L-region,  as 

t 
asserted.    It  is  easily  checked  that  for  each  8  >  0  there  exists  Xo(8)  >  0  such 

that  for  all  X  >  Xo(8)  and  p  ^  e^,  the  5-distance  of  p  to  yi,yj  is  less  than 

X-l-8  (here  we  use  the  fact  that  the  rightmost  point  on  B  is  smooth).  Choose 

8  to  be  the  smallest  B-distance  from  points  in  cq  to  objects  in  R.  Then  the 

above  observation  implies  that  the  B-distance  from  any  point  in  e^  to  Si,  Sj  is 

less  than  its  5-distance  to  any  object  in  R,  as  asserted. 

If  the  rightmost  point  on  B  is  not  smooth  consider  5^  as  defined  above 
for  some  e  >  0.  We  have  the  following  lemma. 

Lemma  3.2:  (a)  For  each  e  >  0  and  p ,  q  i  R^  v^&  have 

ds^(p,q)^  ds(p,q)  ^  (1  +   ^^^^.O,^) 
where  e'  -  0  as  e  -  0  and  A'  >  0  is  the  Euclidean  distance  from  O  to  the 

boundary  of  B. 

(b)  Let  5i  6  L,  ^2  ^  ^  be  two  obstacles.  For  each  e  >  0  let  T^  be  the  locus  of 
points  X  satisfying 

and  similarly  let  F  be  the  locus  of  points  x  satisfying 

ds{x,  ij)  --  J^(.v,  S2). 
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Then  T  is  the  limit  set  (in  the  Hausdorff  topology  of  closed  sets)  of  T^  as 
€^0(Fig3.4).  N  I 

B    CONTOUR 


B, CONTOUR 


'or 


B 


(R) 


-•-Vor^lL) 


Fig.  3,4 

Proof:  (a)  Since  B^  D  B  the  left  inequality  is  obvious.  As  for  the  other 
inequality,  it  is  easy  to  see  that  for  each  e'  >  0  there  exists  8  >  0  such  that 
for  each  e  <  6  we  can  write  the  vector  from  O  to  any  point  on  the  boundary 
of  Z?4  as  X  +  z  where,  x  and  z  are  collinear,  x  ^  Bq  and  llzll  <  e'.  Suppose 
q  ^  p  +  XB^.  Then  we  can  write 

q  —  p  —  \{x  +  z) 
where  x,  z  are  as  above.   But  then  — r-z  i  B  and  we  can  write 


\{x  +  z)  =  X(l  +  ^)x' 

K. 


where 


,   e'     A- 
X+---Z 


'-Y 


e  B 


Hence  d^ip ,  q) 


\(1  +  — )  and  (a)  is  thus  established. 
K 
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(b)  First  we  show  that  if  6„  -  0  and  x„  €  T^^  cxjnverge  to  some  x,  then  x  ^T. 
We  have 

By  (a)  we  then  have,  for  an  appropriate  corresponding  sequence  e'„  -  0, 

€'        -1 

(1  +  -f )     deix^,  s{)  <  ^B  (a:„,  ^O  =  ^^^  (x„,  ^2)  ^  ^fi(^n,  -^2) 

A  «  « 

SO  that  by  continuity  of  cIb  i 

and  by  symmetry  we  obtain  equality,  i.e.  x  ^  F. 

Next  let  ;c  e  r,  that  is  dB{x,s^)  =  Js(x,  ^2),  and  let  8  >  0.  By  properties 
(4),  (9)  of  section  2  it  follows  that  the  disc  Dg  of  radius  5  about  x  must 
contain  two  points  y,  z  such  that 

and 

^fi(2,  ^1)  <  dsiz,  S2)  . 
Hence  if  e  is  sufficiently  small  we  will  also  have 

^fi,0''  -^i)  ~  ^B^y^  ^2)  >  0  >  dej^z,  5i)  -  dsi^z,  S2)  . 
This  and  the  continuity  of  d^^  implies  that  Dg  n  T^  ^  0  and  this  proves  (b). 

Q.E.D. 

Lemma  3.2  now  implies  that  the  contour  C  is  the  limit  set  of  the 
contours  C^,  defined  in  a  similar  manner  for  B^,  as  e  -  0.  This  however 
shows  that  for  every  sufficiently  large  disc  D  about  the  origin 
C^  n  D  ^C  n  D  as  e  -0  and  thus  C  has  the  form  asserted  by  the  lemma. 
Q.E.D. 
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Using  Lemma  3.1  we  can  compute  the  contour  in  the  following  way. 
Following  Kirkpatrick  [Ki]  we  divide  each  cell  of  VorgiL),  VorB{R)  into 
subcells  by  adding  line  segments  which  join  each  vertex  v  of  the  diagram  to 
the  three  points  on  the  three  corresponding  objects  in  L  or  in  R  to  which  v  is 
5-closest.  We  further  refine  the  subcells  by  adding  new  vertices  to  VorB{L) 
(resp.  Vorg{R))  and  joining  these  vertices  by  straight  line  segments  to  their 
fi-closest  points  on  L  (resp.  R)  so  that  each  subcell  will  contain  points  which 
are  5-closest  to  just  one  obstacle  segment  or  corner.  These  line  segments  are 
called  (as  in  [Ki])  spokes.  The  boundary  of  each  resulting  subcell  thus 
consists  of  two  spokes,  an  edge  on  the  5-Voronoi  diagram  (or  possibly  two 
such  edges  if  the  subcell  is  unbounded)  and  a  (possibly  empty)  subsegment  of 
an  obstacle  edge.  Qearly  the  number  of  subcells  and  spokes  is  linear  in  the 
number  of  vertices  of  the  5-Voronoi  diagram  and  the  number  of  obstacle 
corners  (i.e.  0{N  +  n)). 

The  contour  is  now  traced,  starting  at  the  point  ci  defined  in  Lemma  3.1, 
in  the  following  way,  which  resembles  similar  contour  tracing  techniques  as 
described  in  [Ki],  [Sha],  [LS].  Assume  that  we  are  at  a  point  z  on  the  contour 
which  is  an  interior  point  of  a  subcell  U  of  Cs{si)  in  Vorg{L)  and  of  a  subcell 
V  of  Cb{sj)  in  Vorg(R),  and  let  h-,  (resp.  h-^)  be  the  side  (or  comer)  on  Si 
(resp.  Sj)  to  which  z  is  5-closest  (Fig.  3.5(a)).  We  trace  the  contour  portion 
containing  z,  which  is  the  locus  of  points  which  are  equally  fi-close  to  h',  and 
Wj,  until  its  first  intersection  q  with  the  boundary  of  either  U  or  V  (note  that 
the  way  in  which  wc  have  partitioned  the  diagram  into  subcells  and  the 
assumptions  that  we  have  made  on  the  shape  of  B  enable  us  to  find  this 
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(a) 
Fig.  3.5  Contour  tracing  procedure. 


(b) 


intersection  point  in  constant  time).  Suppose  that  the  contour  intersects  the 
boundary  of  V  first  and  enters  into  a  new  subcell  V  of,  say,  Cgis^)  in 
VorsiR)  and  let  wt,  be  the  side  (or  comer)  to  which  all  points  in  V  are  B- 
closest  among  all  objects  in  R.  Then  after  the  intersection  q  the  contour 
continues  along  a  new  subarc  consisting  of  points  equally  5-close  to  w,-  and 
H'jfc.  In  some  cases  (but  only  if  the  leftmost  or  rightmost  point  of  B  is  not 
smooth)  both  subcells  U  and  V  may  be  unbounded  and  the  contour  portion 
being  traced  may  extend  to  infinity  (Fig.  3.5(b)).  In  such  cases  we  know 
from  Lemmas  3.1  and  3.2  that  the  contour  "returns  from  infinity"  from  the 
same  direction  and  we  can  proceed  as  follows.  Suppose  that  the  contour 
portion  being  traced  extends  to  infinity  to  the  left  (resp.  right)  of  the  vertical 
line  separating  L  and  R.  We  then  replace  the  subcell  V  (resp.  V)  by  the 
unbounded  subcell  U'  (resp.  V)  which  is  adjacent  and  lies  clockwise  (resp. 
counterclockwise)  to  it  in  Vorg{L)  (resp.  Vors^R))  and  proceed  as  before  to 
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obtain  the  next  (also  unbounded)  contour  portion. 

We  continue  tracing  the  contour  in  this  manner  until  the  point  C2  defined  in 
Lemma  3.1  is  reached.  In  the  same  manner  we  trace  all  the  other  portions  of 
the  contour.  To  complete  the  construction  of  the  diagram  we  discard  those 
portions  of  the  contour  which  form  the  boundai'y  of  objects  which  were  split 
by  the  vertical  line  separating  L  and  R.  The  complexity  of  this  tracing 
procedure  is  plainly  linear  in  the  number  of  intersections  of  the  contour  with 
edges  and  spokes  of  both  VorB{L),  VorB{R).  The  following  Lemma  proves 
that  the  number  of  such  intersection  points  is  linear  in  the  number  of  objects 
in  L  and  R. 

Lemma  3.3:  (a)  The  contour  may  either  intersect  a  spoke  of  Vorg(L)  or  of 
Vorg{R)  once  or  contain  just  one  subsegment  of  such  a  spoke. 

(b)  The  contour  intersects  edges  of  VorB{L)  or  VorB{R)  at  most  0{m)  times 
where  m  is  the  number  of  objects  in  L  and  R . 

Proof: 

(a)  Follows  from  the  fact  that  each  cell  of  Vorg{S{W))  is  star-shaped. 

(b)  Follows  from  the  fact  that  each  intersection  of  the  contour  with  an  edge 
of  Vor^iL)  or  of  Vorg(R)  is  a  vertex  of  VorB{S(W))  and  there  are  only  0{m) 
such  vertices .    Q.E.D. 

These  observations  yield  the  following  main  theorem. 

Theorem  3.1:  The  5-Voronoi  diagram  for  a  set  S  of  A'  polygonal  convex 
obstacles  consisting  altogether  of  n  sides  can  be  computed  in  time 
0{n  log  A'). 
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Proof:  The  algorithm  consists  of  0{  log  N)  stages.  At  each  stage  we  compute 
fi-Voronoi  diagrams  only  for  windows  which  contain  at  least  one  point  of  P 
and  hence  there  are  only  0{N)  such  "active"  windows  at  each  stage.  The 
number  of  objects  in  each  window  W  (recall  that  these  objects  are  either 
Uw,Dyir  or  intersections  of  objects  in  5  v/ith  W)  is  at  most  two  plus  the 

number  of  points  of  P  within  W.    Summing  over  all  windows  for  which  B- 

i 
Voronoi  diagrams  are  being  computed  in  the  current  stage,  we  see  that  the 

total  number  of  objects  appearing  in  those  windows  is  at  most  the  number  of 
points  in  P  plus  twice  the  number  of  active  windows.  Since  there  are  only 
0{N)  such  active  windows  we  conclude  that  the  number  of  "active"  objects  at 
each  stage  of  the  algorithm  is  0(N).  Therefore  the  total  number  of  wall 
segments  and  corners  in  the  "active"  windows  at  each  stage  is  0(N  +  n)  and 
hence  from  Lemma  3.3  and  the  argument  preceding  it  it  follows  that  the 
computation  done  at  each  stage  is  also  0(N  +  n).  Q.E.D. 
Note:  We  have  assumed  that  the  object  B  is  of  some  simple  form  so  that  we 
can  locate  the  intersections  of  the  contour  with  the  boundary  of  each  subcell 
in  constant  time.  If,  however,  B  is  a  polygonal  object  with  k  sides  then  we 
can  divide  each  cell  into  a  more  refined  set  of  subcells  in  such  a  way  that  the 
5 -Voronoi  edge  on  the  boundary  of  each  subcell  is  a  straight  line  segment. 
The  number  of  subcells  will  then  be  0{kN  +  n)  and  the  algorithm  will  work 
in  time  0{{kN  +  n)  log  A')  (cf.  Property  (6)  of  Section  2). 

4.   Translational  motion  planning  using  Z?-Voronoi  diagrams. 

In   this  final   section   \vc   sliow   how  the  translational  motion  planning 
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problem  described  in  the  introduction  can  be  solved  using  B-Voronoi 
diagrams.  Let  S  =  {si,  .  .  .  ,sj^]  be  a  the  set  of  convex  obstacles  amidst  which 
a  given  convex  body  B  is  allowed  to  translate,  and  let  yorg(5')  be  the  B- 
Voronoi  diagram  of  5.  In  this  section  we  will  impose  the  additional 
requirement  that  the  space  V  bounded  by  the  objects  in  5  is  a  bounded  set. 
This  is  no  real  loss  of  generality  because,  given  S,  it  is  easy  to  see  that  there 
exists  a  large  enough  square  R  (which  can  be  efficiently  computed  from  S) 
such  that  a  free  motion  between  two  positions  p,qiVr\Roib  exist  within 
V  if  and  only  if  it  exists  within  V  C\  R.  Let  us  define,  for  each  point  p  in 
which  (the  reference  point  O  on)  B  can  be  placed,  a  set  r{p)  of  points  on 
Vorg{S)  as  follows.  Suppose  p  ^  V  lies  in  the  interior  of  some  5-cell  Cb{s^. 
Let  X  be  the  B -closest  point  to  p  on  s^.  Then  define  r(p)  to  be  the  intersection 
of  the  ray  from  x  which  passes  through  p  with  the  boundary  of  the  5-cell 
Cb{si).  If  p  6  Vora(5')  then  we  put  r{p)  —  [p].  By  the  boundedness  of  V  and 
by  property  (2)  of  section  2,  rip)  is  always  nonempty,  and  is  either  a  single 
point  or  a  straight  line  segment  (see  Fig  4.1(a)).  Let  VorB'{S)  be  a 
subdiagram  of  Vorg{S)  which  is  composed  of  all  points  of  Vorg{S)  whose  B- 
closeness  to  S  is  greater  than  one.  We  can  easily  compute  Vorg'{S)  from 
'Vorg{S)  in  linear  time  by  considering  each  subccll  of  Vorg{S)  separately.  We 
have  the  following  theorem. 

Theorem  4.1:  Let  p ,  qho.  two  positions  of  B  with  the  same  fixed  orientation. 
Then  there  exists  a  continuous  obstacle-avoiding  purely  translational  motion 
of  B  from  p  to  ^  if  and  only  if  r{p)  and  r{q)  belong  to  the  same  connected 
component  of  'Vorg'{S). 
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Proof:  Suppose  first  that  r{p)  and  r{q)  lie  in  the  same  connected  component 
of  Vorg^S).  Then  the  required  motion  of  B  from  p  to  ^  is  obtained  by  first 
moving  (the  reference  point  O  on)  B  along  the  straight  segment  from  p  to  a 
point  in  r{p),  then  to  a  point  in  r{q)  along  a  path  in  Vor^'iS)  which  connects 
these  two  sets,  and  finally  along  the  straight  segment  from  r{q)  back  to  q. 
(Note  that  the  initial  and  final  portions  of  that  motion  are  collision  free 
because  the  5-distance  to  5  keeps  increasing  as  we  move  from  p  (resp.  q) 
along  the  straight  segment  to  r{p)  (resp.  r{q)).)  Conversely,  let 
v(r),  0  <  /  <1,  be  a  continuous  path  in  V  that  O  traces  as  B  is  translated 
from  p  -  v(0)  to  ^  =  v(l)  without  colliding  with  any  obstacle.  It  is  clear 
that  dB{v{t),S)  >  1  for  every  t  since  otherwise  B  would  intersect  an  obstacle 
during  that  motion.  By  definition  of  the  5-Voronoi  diagram  we  have  for 
each  /,  di;{r{v{t)),S)  >  ds{v{t),S).  Thus  the  union  K  of  all  the  sets 
r(v(r)),  0  <  r  <1,  is  contained  in  Vorg'iS).  We  claim  that  K  is  locally 
connected  at  a  sufficiently  small  neighborhood  of  each  r{v{t))  for 
v(r)  i  Vorg{S).  We  can  also  assume  without  loss  of  generality  that 
v(r)  e  Vorg{S)  only  finitely  many  times,  and  we  will  show  how  to  complete  K 
for  each  such  t  to  make  it  fully  connected  while  still  being  contained  in 
Vors'(S). 

To  show  the  above  claim  it  suffices  to  show  that  if  p„  -p,  with  p  i  Vorg{S), 
then  there  exists  a  sequence  h'„  i  r(p„)  which  converges  to  some  point  in 
r(p)  (Fig.  4.1(a)).  Let  p^-p  and  for  each  n  >  1  let  h'„  6  r{p„)  be  an 
arbitrary  point  such  that  the  sequence  u'„  converges  to  some  point  w.  By 
passing  to  subsequences  if  necessary,  we  may  assume  that  there  exist  two 
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objects  si,  S2  in  S  such  that  for  each  n  >  1,  p,p^  ^  Cb{si), 
^ai^n^si)  =  dB{wn,S2)  ^  ^bCw^,-^)  for  J  i^  si,S2,  the  straight  segments  j:„p„w„ 
are  all  contained  in  Csisi),  where  x„  is  the  point  on  s-^  which  is  B-closest  to 
p„  (hence  also  to  >v„;  cf.  Property  (2)  of  Section  2),  and  x„  converge  to  some 
X  €  5i.  By  continuity  these  properties  will  also  hold  for  x,  p,  and  >v.  But 
then  X  must  be  the  B-closest  point  to  p  and  to  w  on  Si,  w  must  belong  to 
VorB{S),  and  the  intersection  of  the  ray  from  x  through  p  with  the  boundary 
of  Cb{si)  must  contain  w.  By  definition  w  C  r(p)  and  the  above  claim  is  thus 
established.  If  p  ^  Vorg(5)  (Fig.  4.1(b)),  then  the  above  arguments  remain 
the  same,  except  that  w  may  not  belong  to  r(p).  However,  since  both  p  and 
w  belong  to  Cg(5'i)  and  p  lies  on  the  segment  from  x  to  w,  it  follows  from 
Property  (2)  of  Section  2  that  the  segment  pw  is  contained  in  Vor^iS) 
(actually  in  Vorg'^S)  as  is  easily  checked).  Hence  if  we  add  to  A'  a  finite 
number  of  such  straight  segments  for  points  at  which  v(/)  intersects  Vorg(S), 
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we  obtain  a  connected  subset  of  Vorg'^S)  containing  both  r(p)  and  r{q).  This 
completes  the  proof  of  the  Theorem.   Q.E.D. 

Remarks:  (1)  Note  that  Vorg'{S)  may  have  vertices  with  only  one  incident 
edge;  these  are  points  p  whose  B-distance  to  S  is  exactly  1.  In  the 
terminology    of    [KS],    these    points    are    precisely    the    points    of    local 

nonconvexity  on  the  boundary  of  K  =   U  (si  —  B).  Hence  the  5-Voronoi 

diagram  provides  us  with  an  alternative  (and  more  efficient)  method  of 
calculating  these  locally  nonconvex  comers  of  K,  and  thus  leads  to  a  more 
efficient  calculation  of  K  itself  (see  [KS],  [KLPS]  for  more  detail). 

(2)  Another  useful  property  of  VorB{S)  is  that  it  can  be  used  for  planning 
translational  collision-free  motion  of  any  scalar  multiple  \B  of  B.  For  each 
given  X  we  only  need  to  perform  a  linear  time  pre-processing  of  Vorg{S)  to 
produce  a  truncated  subdiagram  similar  to  Vorg'{S)  and  then  plan  each 
required  motion  by  an  appropriate  (linear-time)  graph  searching  through  that 
subdiagram  (cf.  [OY]  for  a  similar  remark  concerning  motion  planning  for 
circles  of  arbitrary  size). 
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